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ABSTRACT 

We present observational data for two main components (S255IR and S255N) of the S255 high mass 
star forming region in continuum and molecular lines obtained at 1.3 mm and 1.1 mm with the SMA, 
at 1.3 cm with the VLA and at 23 and 50 cm with the GMRT. The angular resolution was from ~ 2" to 
~ 5" for all instruments. With the SMA we detected a total of about 50 spectral lines of 20 different 
molecules (including isotopologues). About half of the lines and half of the species (in particular 
N 2 H+, SiO, C 34 S, DCN, DNC, DCO+, HC 3 N, H 2 CO, H 2 CS, S0 2 ) have not been previously reported 
in S255IR and partly in S255N at high angular resolution. Our data reveal several new clumps in 
the S255IR and S255N areas by their millimeter wave continuum emission. Masses of these clumps 
are estimated at a few solar masses. The line widths greatly exceed expected thermal widths. These 
clumps have practically no association with NIR or radio continuum sources, implying a very early 
stage of evolution. At the same time, our SiO data indicate the presence of high-velocity outflows 
related to some of these clumps. In some cases, strong molecular emission at velocities of the quiescent 
gas has no detectable counterpart in the continuum. We discuss the main features of the distribution 
of NH3, N 2 H + , and deuterated molecules. We estimate properties of decimeter wave radio continuum 
sources and their relationship with the molecular material. 

Subject headings: astrochemistry - HII regions - instrumentation: interferometers - ISM: clouds - 
ISM: molecules - stars: formation 



1. INTRODUCTION 

Despite their great importance in almost all areas of 
astronomy, the formation of stars larger than 8-10 M 
is still poorly understood. In part this is because high 
mass star formation (HMSF) regions are more distant, 
more active, and shorter lived than their low-mass coun- 
terparts. Nevertheless, an enhanced attention has been 
paid recently to the earliest phases of massive star forma- 
tion. So far, few young HMSF cores have been studied 
in detail and it is important to extend the list of such ob- 
jects, especially considering the large variety of processes 
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that probably lead to the formation of a massive star. 

S255 is an H 11 region associated with a dense core at 
a commonly accepted d istance of 2.5 kpc (Russeil et al. 
2007} IQjha et al.| | 20 11 [) which we ad opt here too. We 
note, however, that Rygl et aT] ( |2010[ ) report a distance 
of 1.6 kpc based on trigonometric parallax measurements 
of methanol masers. The core consists of two main com- 
ponents (S255IR and S255N) separated by slightly over 
V . Our team has previously acquired molecular line 
observations using single-dish instruments (OSO-20m, 
IRAM-30m, NRA O-12m). With angular r esolution from 
about V to 10" (Zinchenko et al. 2009), we obtained 
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AV ~ 2 kms -1 for both components. While both com- 
ponents show evidence for cool, massive clumps, their 
evolutionary states appear to be quite different. S255IR 
is bright (> 70 Jy) in all the Midcourse Space Exper- 
iment (MSX) ban ds and contains a near-IR cluster o f 
early-B-type stars ( |Howard et al.|1997lfltoh et alJ200l| ), 
a cluster of compact Hii regions ( pneil & Bally 1986), 
and a wea lth of complex H2 emission features ( |lVliralles| 
et al.||1997D . In contrast, S255N (also called Sh2-255 
FIETand G192.6 0-MM1) contains a single cometary UC 
Hii region (e.g., [Kurtz et al.|1994| ), and is u ndetected by 
MSX at w avelengths shorter than 21 /am ( jCrowther fc| 
Conti|[2QQ3 ). VLA and SMA observations in the contin- 



uum, several molecular lines and water maser emission 
( Cyganowski et al.|2QQ7 ) indicate the presence of a mas- 
sive protocluster in this region. The chemical composi- 
tion of S255IR and S255N also appear sig nificantly dif- 
ferent flLintott et al.|2005| [Zinchenko et al. 12009] ). While 
the CS and H(JJN abundances are very similar, the abun- 
dances of NH3, N2H + , HCO + and some other molecules 
in these components differ significantly. 

Recently, both components were studied with the SMA 
at 1.3 mm in CO 13 CO, C 18 Q CH 3OH, CH 3 CN and 
some other lines ( |Wang et al. 2011). These observa- 
tions revealed 3 continuum clumps m the S255IR area 
and high-velocity collimated outflows in both regions. 
The star p opulation in this complex was studied most 
recently by |Ojha et al.| ( 2011 ) on the basis of optical and 
NIR observations. They found a number of new YSO 
candidates with a large spread in ages, indicating a sce- 
nario of induced star formation. 

In general, this star forming complex represents an ex- 
cellent laboratory for studies of different stages of mas- 
sive star formation. Our goal is to investigate further the 
structure, physical properties and chemistry of this area 
on small scales by observations of important molecular 
tracers like NH3 (which is a convenient "thermometer" 
for dense molecular cores and provides i nformation on 
their density and kinematics too 
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|l983||Walmsley fc Ungerechts|l983[ ), N 2 H+, SiO, deuter- 
ated molecules, etc. In addition we want to investigate in 
more detail the distribution and properties of ionized gas 
in this region in order to trace the interaction of massive 
YSOs and outflows with the surrounding medium. 

In this paper we present observational data obtained 
at 1.3 mm and 1.1 mm with the SMA in the compact 
configuration (both in various lines and in continuum), 
with the VLA at 1.3 cm (in ammonia lines) and with the 
GMRT at 23 and 50 cm in continuum. We also discuss 
the general structure and kinematics of the complex and 
derive basic physical parameters of the observed features. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. SMA 

The S255 complex, including S255 N and S255 IR, 
was observed with the Submillimeter Array (SMA) 
in its compact configuration on 2010 April 7th 
and April 14th at 224.8 GHz and 283.9 GHz, re- 
spectively. Two fields, one centered at S255 N 
(06H2 m 53 s .669, 18°00'26".903) and one centered at 
S255 IR (06H2 m 53 s .800, 17 o 59'22".097) were inter- 
leaved throughout the observations. The primary HPBW 
of the SMA antennas is 55" — 45" at these frequencies. 



Typical system temperatures on source were between 
90 K and 160 K for the 224.8 GHz track, and between 
100 K and 200 K for the 283.9 GHz track. The result- 
ing ^v-coverage ranges from 6 kA to 50 kA and 8 kA 
to 65 kA for the 224.8 GHz and 283.9 GHz tracks, re- 
spectively. 3C273 was used as the bandpass calibrator 
for both tracks, and 0530+135 as well as 0750+125 were 
used as the complex phase and amplitude gain calibrator. 
The gain calibrator flux scale, calibrated against Mars, 
was found to be consistent within 5% with the SMA cal- 
ibrator database and estimated to be accurate within 
20%. In the 224.8 GHz observation, a total of 8 GHz 
(216.8-220.8 GHz in the LSB and 228.8-232.8 GHz in 
the USB) was observed with the SMA bandwidth dou- 
bling correlator configuration. A number of different 
spectral resolutions, with a maximum resolution of 0.53 
km s -1 , were employed for different molecular lines. The 
283.9 GHz observation covered 277.9-279.9 GHz in the 
LSB and 287.9-289.9 GHz in the USB. The data calibra- 



tion was carried out with the IDL superset MIR (Scov- 
ille et al. 1993), and s ubsequent imagin g and analysis 
were done m MIRIAD flSault et al ; l|1995| >. With robust 
weighting for the continuum and line data, the synthe- 
sized beam sizes are approximately 3'/ 8 x 3'/0 (PA 60°) at 
1.3 mm and 2 / /9x2 / /6 (PA 40°) at 1.1 mm. The noise level 
is ~ 8 mJy beam -1 in the 1.1 mm continuum image and 
roughly 60 mJybeam -1 and 100 mJybeam -1 in the line 
cubes at 2.0 kms -1 spectral resolution at 220-230 GHz 
and 280-290 GHz, respectively. 

2.2. VLA 

The VLA observations (program AZ186) were made 
on 2009 November 13 (S255N) and 21 (S255IR) with 
the array in the D-configuration. The pointing cen- 
ters for S255N and S255IR in J2000 coordinates were 
06H2 m 53 s .7 +18°00 / 27 // and 06 /l 12 m 54 s .O +17°59 / 22 // , 
respectively. At the observing frequency the VLA pri- 
mary beam is approximately 2 arcmin. The correla- 
tor was in the 1IF mode, measuring right circular po- 
larization with a 6.25 MHz bandwidth and 127 chan- 
nels of 48.8 kHz each. The ammonia (1,1) and (2,2) 
lines were observed, with rest frequencies of 23,693.7955 
and 23,721.9336 MHz. 3C147 was observed as the flux 
calibrator, and J0530+135 was observed as the phase 
and bandpass calibrator. The assumed flux density for 
3C147 was 1.700 Jy and the bootstrapped flux density for 
J0530+135 was 1.439 + 0.015 Jy. The angular resolution 
of the image cubes was 2 / /61 x 2 / /51 (PA -65.6°). The 
approximate rms noise of the cube was 3.5 mJy beam -1 . 
The on-source time was approximately one half hour in 
each transition. 

2.3. GMRT 

The radio continuum interferometric observations at 
610 MHz and 1280 MHz were carried out on 2009 May 
08 and 2009 September 08, using the Giant Metrewave 
Radio Telescope (GMRT) array. The GMRT has a "Y"- 
shaped hybrid configuration of 30 antennas, each 45 m in 
diameter. There are six antennas along each of the three 
arms (with arm length of ~ 14 km). The remaining 
twelve antennas are located in a random and compact 
(~ 1km x 1km) arrangement near the center. Details 
of the GMRT antennas and their configurations can be 
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found in Swamp et al. (1991). The largest angular scales 
to which the GMKT is sensitive are 8' and 17' at 1280 
and 610 MHz, respectively. 

Flux and phase calibrators were observed to derive the 
phase and amplitude gains of the antennas. 3C147 was 
used as the flux calibrator, while the phase calibrators 
were 0503+020 and 0632+103 at 610 MHz and 1280 
MHz, respectively. The data analysis was done using 
AIPS. Corrupted data were identified and flagged. The 
estimated uncertainty of the flux calibration is within 
7% at both frequencies. Images of the field were formed 
by the Fourier inversion and cleaning algorithm task 
IMAGR. Several iterations of self-calibration were car- 
ried out to remove the residual effects of atmospheric and 
ionospheric phase corruption and to improve the quality 
of the maps. The high resolution images (6'/ 4 x 4'/ 8 at 
610 MHz and 5'/2 x 3'/8 at 1280 MHz) were made with a 
Briggs weighting function halfway between uniform and 
natural weighting (robust factor = -0.5 to 0), which is a 
good compromise between angular resolution and sensi- 
tivity. 

3. RESULTS 

With the SMA we detected about 50 spectral lines of 
20 different molecules (including isotopologues). About 
half of the lines and half of the species (in particular 
12 CO, 13 CO, C 18 0, CH 3 OH, CH 3 CN, 13 CS, HNCO, SO, 
P CS) were previously o bserved in this area with the SMA 
by I Wang et al.| (2011). We do not discuss these lines 
in detail; rather, we compare their results with o ur line 
detections. Our additional to Wang et al. ( 2011| ) obser- 
vations include lines of N 2 H+, SiO, C^ 4 S, DCN, DNC, 
DCO+, HC 3 N, H 2 CO, H 2 CS, S0 2 as well as additional 
transitions of CH 3 OH, HNCO and OCS. A list of these 
lines, including their frequencies and lower energy levels, 
is given in Table fl] Some of them (marked by asterisk 
m Table [l]) were observed earlier in S255N at the SMA 
by Cyganowski et al. (2007). Nevertheless, they have 
not been observed in S2551R and our data show many 
new features in these lines in S255N, which are presented 
and discussed below. In addition, we observed the (1,1) 
and (2,2) ammonia lines with the VLA. Here we particu- 
larly focus on several key species, including N 2 H + , NH3, 
deuterated molecules, SiO, methanol, etc. It is worth 
noting that we did not detect some important lines in 
our band, in particular, N 2 D + and radio recombination 
lines. 

We present the results in the form of maps as well as 
spectra and line parameters at selected positions. For 
continuum observations we give positions, flux densities, 
and size estimates of the continuum sources. In general 
the results are presented separately for the S255IR and 
S255N areas. 

3.1. Millimeter wave continuum 

Maps of the S255IR and S255N areas in 1.1 mm contin- 
uum overlaid on NIR K-band images flOjha et al.|[20TT| ) 
are shown in Fig. [I] The two brightest continuum sources 
visible in the S255IR map have been reported and dis- 
cussed by [Wang et al] ffim\ (S255IR-SMA1 and S255IR- 
SMA2 in their nomenclature) . In addition, they identi- 
fied the S255IR-SMA3 clump which is indistinguishable 
from S255IR-SMA1 in our data (due to our lower angu- 
lar resolution). On the other hand, our map shows an 



Table 1 

List of molecular transitions observed at the SMA in 
either S255IR or S255N in addition to those observed by 
|Wang et al~1(|2011|>. 



Molecule 


Transition 


Frequency 
(GHz) 


Ei 
(cm- 1 ) 




CH3OH 


5i - 4 2 E 


216.945521 


31.596 






61 - 7 2 A- 


217.299205 


252.644 






4 2 -3i E 


218.440050 


24.310 


* 




10 2 - 9 3 A+ 


232.418571 


107.208 






9_i - 8 E 


278.304575 


67.150 






14 4 - 15 3 E 


278.599037 


226.775 






11 2 - 10 3 A- 


279.351887 


123.338 






60 - 5 E 


289.939386 


33.316 




SiO 


5-4 


217.104984 


14.484 


* 


c 34 s 


6-5 


289.209068 


24.119 




DCN 


3-2 


217.238530 


7.246 






4-3 


289.644907 


14.493 




DNC 


3-2 


228.910471 


7.636 




DCO+ 


4-3 


288.143858 


14.418 




HNCO 


10i,io -9i, 9 


218.981170 


62.949 




HC 3 N 


24-23 


218.324711 


83.755 


* 


H 2 CS 


81,7 — 7i,6 


278.887680 


41.733 




H 2 CO 


3o,3 — 2o, 2 


218.222195 


7.286 


* 




3 2 , 2 — 22,1 


218.475642 


40.040 


* 




32,1 — 22,0 


218.760071 


40.043 


* 


N 2 H+ 


3-2 


279.511760 


9.324 




OCS 


18-17 


218.903357 


62.070 






23-22 


279.685318 


102.634 




S0 2 


18l,17 — 172,16 


288.519996 


103.712 





Observed earlier at the SMA in S255N by Cyganowski 
[eTaL] (2007) 

additional, very bright clump about 10" further south. 
We designate this clump as S255IR-SMA4. 
In the S255N area, in addition to the continuum 



clumps rep orted earlier (Cyganowski et al. 2007 [Wang 
et al. | 2011 ) our 1.1 mm continuum map shows three new 



features that are also detected in molecular emission (see 
below). We designate these as S255N-SMA4, SMA5 and 
SMA6; their parameters are given in Table[| 

3.2. Molecular emission 

Our data show that the intensity distributions for 
many molecular species are quite different from the con- 
tinuum maps. In many cases, the molecular emission 
peaks do not coincide with the continuum peaks and 
there is also molecular emission without corresponding 
continuum emission. Below we describe in more detail 
the emission features of several key species. 

3.2.1. NH 3 

The ammonia emission i n the S255IR area is r ather 
weak (in accordance with |Zinchenko et al." 1997) and 
shows significantly different distributions in the NH3 
(1,1) and (2,2) transitions (Fig. [2]). The (1,1) emission 
peaks at the SMA2 clump while the (2,2) emission peak 
coincides with SMA1. This indicates that the temper- 
ature of the SM Al clump is much h igher than that of 
SMA2. Indeed, J Wang et aT7] fl201ip derived a kinetic 
temperature of 150 K tor SMAlTloased on CH3CN data, 
implying the presence of a hot core. The main emission 
peak towards SMAl is observed at Vlsr ~ 4.8 kms -1 . 
There is a hint on another peak at Vlsr ~ 9 kms -1 . 

The ammonia emission is relatively strong in S255N. 
In both the (1,1) and (2,2) transitions it extends far to 
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Figure 1. Maps of the 1.1 mm continuum e mission in the S255 IR (left panel) and S255N (right panel) areas (contours) overlaid on the 
K band NIR images (in logarithmic scale) from |Ojha et al.| ( |201l| . The contour levels are (2, 2.5, 3, 3.5, 4, 5, 6, 7, 8, 10)x25 mJybeam -1 
for S255IR and (-4, -3.5, -3, -2.5, -2, 2, 2.5, 3, 3.5, 4, 5, 6, 7, 8, 10, 12, 15, 20, 25) x 12.5 mJybeam- 1 for S255N. The dashed contours 
show negative features due to the missing flux. The SMA beam is shown in the lower left corner of both panels. 



Table 2 

Positions, flux densities, deconvolved angular sizes and position angles of the 
millimeter wave continuum sources at 284 GHz as well as flux densities at 
225 GHz (the last column) measured with the SMA. 



Name 


a(2000) 


5(2000) 


#284 


$max 


$min 


P.A. 


#225 




(h m s) 


(° ' ") 


(Jy) 


II 


II 


(°) 


(Jy) 


S255IR-SMA1 


6:12:53.98 


17:59:23.1 


0.45 


1.5 


0.7 


-14 


0.29 


S255IR-SMA2 


6:12:53.75 


17:59:25.7 


0.43 


2.9 


2.2 


11 


0.30 


S255IR-SMA4 


6:12:54.03 


17:59:11.5 


0.31 


4.4 


3.5 


-50 


0.08 


S255N-SMA1 


6:12:53.64 


18:00:26.8 


0.75 


3.0 


2.0 


4 


0.54 


S255N-SMA2 


6:12:52.95 


18:00:31.7 


0.05 


2.5 


2.1 


-10 


0.06 


S255N-SMA3 


6:12:53.62 


18:00:18.7 


0.05 


3.3 


2.9 


-31 


< 0.03 


S255N-SMA4 


6:12:54.32 


18:00:41.0 


0.16 


3.8 


0.7 


-12 


0.04 


S255N-SMA5 


6:12:54.00 


18:00:13.2 


0.15 


4.5 


3.5 


49 


0.04 


S255N-SMA6 


6:12:53.04 


18:00:03.4 


0.14 


7.9 


1.3 


-40 


< 0.03 



i i 1 i i i i 1 i i 
NH 3 (1,1) 


i i i i i i 
9 S255IR 


i i i i i i i i i i 
NH 3 (2,2) 


i i i i i i 
S255IR 


N 2 H + (3-Gp S255IR 
() 




O 

o ) 

o 

- [4.5,10.7] km s" 1 

i i i i - 1 i i i i i 




.i.i.i 


- [4.5,10.7] km s" 1 

. . i . . . . i i i 




. 1 . , . . 


^0)N 2 H + (1) 
N 2 H + (2)^P^ 

- km s~ 1 

. t i t . i i i i i i t i . i I" d i 


5 h 12 m 55 s 54 s 


53 s 


54 s 


53 s 


54 s 53 s 



a (2000) 

Figure 2. Maps of the NH3 (1,1), NH3 (2,2) and N2H+ (3-2) integrated line emission in the S255IR area (contours) overlaid on the 1.1 mm 
continuum image. The velocity ranges are indicated at the bottom of the plots. The contour levels are (-3, 3, 5, 7)xl0 mJybeam -1 kms -1 
for NH3 and (-3, 3, 5, 7, 9, 11) x 1.4 Jybeam -1 kms -1 for N2H+. The dashed contours show negative features due to the missing flux. 
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Figure 3. Maps of the NH 3 (1,1), NH 3 (2,2) and N 2 H+ (3-2) integrated line emission in the S255N area (contours) overlaid on the 1.1 mm 
continuum image. The velocity ranges are indicated at the bottom of the plots. The contour levels are (-3, 3, 5, 7)xl0 mJybeam -1 kms -1 
for NH3 and (-3, 3, 5, 7)xl.4 Jy beam -1 kms -1 for N2H+. The dashed contours show negative features due to the missing flux. The 
dash-dotted (red in the color version) and dotted (blue in the color version) contours show the red-shifted and blue-shifted high-velocity 
CO emission, respectively. 
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V LS r (km/s) 

Figure 4. The spectra of the NH3 (1,1) and (2,2) emission from 
the S255N-NH3 area averaged over the region 2 / /5 x 3'/ 5 centered at 
R.A.(2000) = 06^12 m 53 s .39, Dec.(2000) = 18°00 / 41 // .0 (in units of 
brightness temperature — the conversion factor is approximately 
320 K per Jy beam" 1 ). The (2,2) spectrum is shifted by -4 K 
along the vertical axis for clarity. 

the north from the main continuum clump (Fig. [3]). This 
extension has no counterpart in the continuum (at any 
observed wavelength) and almost no counterpart in other 
molecular lines (only N2H + and perhaps DCO + show 
similar features, although less extended — see below). 
The (1,1) and (2,2) distributions in this northern exten- 
sion are somewhat different; we discuss this in Section [4j 
We shall designate this northern source of ammonia emis- 
sion as S255N-NH 3 . The (1,1) and (2,2) spectra from this 
area are presented in Fig. [4 They clearly show two nar- 
row (linewidths of about lkrns" 1 ) velocity components 
(at approximately 10.1 and 8.2 kms -1 ). The two ve- 
locity components have significantly different line ratios 
[/(2,2)//(l,l) - 0.2 and - 0.4 for the 8 and 10 kms" 1 
components, respectively], implying different tempera- 
tures for these two kinematic components. The temper- 
ature estimates are given in Section |4j 

Another peak of ammonia emission — the strongest 
one in the (2,2) transition — is displaced by a few arc sec- 
onds to the south from the main continuum source. Only 
the 10 kms -1 kinematic component is readily seen here, 
although there is a hint on the presence of the 8 kms -1 
component in the (2,2) line (at about 1.5<r level). If real, 
it would imply a significant increase of the temperature 
of this component in this area since the 7(2, 2)/I(l, 1) 
ratio appears to be > 1. 

3.2.2. N 2 H+ 

N2H + J = 3 — 2 emission is quite strong in both S255IR 
and S255 N. The integrated intensity maps are shown in 
Figs. |2|3j To better characterize the struct ure and kine- 
matics of these regions, we present in Figs. 5|6 the chan- 
nel maps of this line with 2 kms -1 resolution. 

In S255IR, while there is no detectable N2H + emission 
directly associated with the SMA1 clump, several N2H + 
emission peaks can be seen in this general area (Fig. 
Two of these are related to the SMA2 and SMA4 clumps, 
although the N2H + peaks do not coincide exactly with 
the continuum peaks. Two other peaks have no associ- 
ated continuum. One of these is seen as an extension to 
the west from the SMA2 clump and is apparently related 



to the CH3OH emission feature noticed by |Wang et al. 
(2011). This feature is seen in several other methanol 
lines in our data. NH 3 (1,1) emission is also observed 
here, as are some other lines (e.g., H2CO). We shall des- 
ignate this source as S255IR-N2H + (1). Another quite 
strong N2H + peak, shifted to the SE from SMA1, has no 
counterpart in the continuum or in any observed lines. 
For this emission feature we shall use the name S255IR- 
N2H + (2). Both features are marked in Fig. [2] 

In the S255N area the N2H + emission shows an exten- 
sion in the same direction as NH3 (Fig. [3|. The main 
peak of the N2H + emission is displaced from the contin- 
uum peak (similar to the ammonia) and has an exten- 
sion roughly in the direction of the outflow described by 
Wang et al. (2011). There is one more N2H + emission 
peak about e 2F r io the SSW from S255N-SMA1; it is as- 
sociated with weak continuum emission (S255N-SMA6) 
and is traced by several other molecular lines (see be- 
low). The spectra of N2H + and other molecules from 
S255N-SMA1 are shown in Fig. [f| 

3.2.3. DCN, DNC and DCO + 

Our data set contains lines of several deuterated 
molecules: DCN J = 3-2 and J = 4-3, DNC J = 3-2, 
and DCO + J = 4 — 3. DCN emission is quite strong 
in both S255IR and S255N. In S255IR it is observed in 
both main clumps, SMA1 and SMA2, although its peak 
is shifted from the continuum peak in SMA1 (Fig. [8]). 
The DCN J = 4 - 3 line is also seen in SMA4. In S255N 
the DCN emission coincides with the SMA1 continuum 
clump, having an extension to the north (Fig. [9]). 

The DNC emission is significantly weaker than DCN. 
We did not detect DNC in the S255IR area. In S255N 
it is detected, but the peak is shifted from the DCN and 
continuum peaks. 

DCO + emission is seen in both areas. Its distribution 
is significantly different from that of DCN, and quali- 
tatively resembles that of N2H + . In particular, both 
species avoid ionized regions, and both are observed to- 
wards the S255N SMA6 clump (in contrast to most other 
molecules) . 

3.2.4. SiO 

As expected, the bulk of the SiO emission in both 
S255IR and S255N is apparently associated with out- 
flows. In S255IR the strongest SiO emission is observed 
towards the SMA1 clump (Fig. 10 ) which is probably 
the driving source of the outflow ( |Wang et al. 2011). 
The peak is shifted from the continuum peak in the di- 
rection of the red-shifted outflow lobe. The central ve- 
locity is significantly higher than for most of the other 
lines (Fig. [7|). There is also blue-shifted high- velocity SiO 
emission. K discussion of the high-velocity emission in 
various lines is given in Section [4j 

S255IR also shows narrow-line SiO emission near 
SMA2 and SMA4. The possible association of this emis- 
sion with other features is discussed b elow. T he details 
of SiO kinematics can be seen in Figs. 12|13 In partic- 
ular we see SiO emission from the blue-shifted lobe of 



the extended outflow described by |Wang et al. (2011). 
However, there is no SiO emission from the red- shitted 
lobe of this outflow. 

In S255N (Fig. [TT| ) the main peak of the SiO emission 
coincides with theHSW (blue-shifted) lobe of the molec- 
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ular outflow traced in CO ( Wang et al.||2011 ) and other 
lines. The SiO profile shows an extended blue wing. 
In addition, there are very broad SiO emission features 
clearly seen in the SiO channel maps (Fig. [l3]) associ- 
ated with the SMA3 and SMA5 clumps. In the vicin- 
ity of the SMA3 clump, the SiO emission spans from 
about -45 kins" 1 to about +65 kms -1 (Fig. 14). The 
blue-shifted and red-shifted parts of the emission are spa- 
tially separated as can be easily seen in the channel maps 
(Fig. 13) and more clearly in Fig. 15 All this indicates 
the presence of a high-velocity outnow associated with 
this clump. This outflow is also traced in high- velocity 
CO e mission (Fig. [3] here and Fig. 10 in Wang et al. 
2011] ). 

The total width of the SiO emission associated with 
the SMA5 clump is about 50 kms" 1 (Fig. [Tsb _ The 



blue-shifted and red-shifted SiO emissions are spatially 
separated (Fig. [l7|) . This suggests that another outflow 
is related to this clump. The driving source of the out- 
flow is shifted by a few arcseconds to the south from the 
SMA5 continuum peak. 

3.2.5. CH 3 OH 

In our band we observed several methanol lines of dif- 
ferent excitation. Maps of the methanol emission in the 
lines detected in the S255IR and S255N areas are pre- 
sented in Figs. |18|19 



The main features of the methanol emission in the cen- 
tra l parts of S255IR an d S255N are the same as reported 
by jWang et al.| p011). The highest excitation lines are 
detected only in I3255IR-SMA1. Several methanol lines 
trace some of our new continuum clumps, in parti cular 
S255IR-SMA4, S255N-SMA4 and S255N-SMA5. |Wanj| 



et al. (|201l| reported an emission in one of the methanol 
from 



lines from the area which we designate here as S255IR- 
N 2 H + (1). Our data show an emission here in 3 other 
methanol lines. 

Methanol lines from S255IR-SMA1 and S255N-SMA1 
show wings, although less extended than in SiO. We do 
not see high velocity methanol emission from S255N- 
SMA3 and SMA5 where we report SiO outflows. 

An inspection of the presented maps indicates apparent 
peculiarities in the methanol excitation. For example 
the low excitation 5i — 42 E line is much stronger in 
S255-SMA1 than in SMA2, in contrast to other lines of 
comparable excitation. The same can be said about the 
9_i — 8o E line in the S255N area (the component shifted 
by a few arcseconds to north-east from S255N-SMA1). 
These a nd some other methanol lines in our band can be 
masme; ( |Sobolev|1993l|Cragg et al.|2005| [Voronkov et al.| 
2012 and A. Sobolev, private communication) . 



3.2.6. Other lines 

Our data set contains several other important molecu- 
lar lines including C 34 S, HC 3 N, SO, S0 2 , H 2 CO, H 2 CS, 
CH 3 CN, OCS, and HNCO. Most of them have not been 
observed with high angular resolution in both these ob- 
jects before now. Some of these lines, in particular those 
with relatively high excitation energy (e.g., HC3N and 
SO2), are detected almost exclusively towards the hot 
core in S255IR-SMA1. Other molec ules, (e.g., SO) trace 
the outflows. In Figs. 



It traces most of the new continuum clumps in both 
S255IR and S255N, and also several other features, in- 
cluding S255IR-N 2 H+(1). 

It is worth noting that our C 34 S J = 6 — 5 map of the 
S255IR area looks very different from t he 13 CS J = 5 — 4 
map presented by |Wang et al.| (2011). Our data show 
the C 34 S emission is of similar intensity in both S255IR- 
SMA1 and S255IR-SMA2, and approxi matel y coincident 
with the d ust continuum peak s (Figs. 7|10 ). The 13 CS 
emission of |Wang et al.| ( |2011| ) is seen only in S255IR- 
SMA1 and is significantly shifted from the continuum 
peak. Interestingly, our 13 CS observations of the same 
line show the intensity distribution very similar to C 34 S 
J = 6 — 5. The only reasonable explanation is that the 
bulk of the 13 CS emission i s resolved out by the smaller 
beam of | Wang et al.| ( |2011[ ) (which is about 2 times nar- 
rower than our beam) and they see only a compact struc- 
ture at the eastern edge of S255IR-SMA1. 

3.3. Decimeter and centimeter wave radio continuum 

emission 



10|11|20|21| we present maps of 
the emission of C 34 S, H0 3 N, H 2 CO, H 2 CS, and S0 2 . In 
particular, fairly strong emission is observed in H 2 CO. 



In Fig. 22 we show a large-scale 610 MHz map of the 
region, obtained with the GMRT. It covers the evolved 
H11 regions S255, S256 and S257. Two compact sources 
in the center are associated with the S255IR and S255N 
star-forming cores. 

In Fig. [23] we present 1280 MHz GMRT continuum 
maps of S255IR and S255N. A comparison of the original 
GMRT maps with high-frequency VLA archival maps (in 
particular, at 15 GHz) shows a positional shift of a few 
arc seconds between the maps. The reasons for this shift 
are unclear but probably are related to phase calibration 
at the GMRT. We made positional corrections to the 
GMRT maps to compensate for the shifts. Thus, the 
absolute coordinates in these maps are not reliable to 
more than a few arc seconds. 

The main components seen in th e GMRT maps were 
identified by Snell & Bally ( 1986| ) and we follow their 
nomenclature here tor the radio continuum sources (in 
particular S255-1 corresponds to S255N). Our data ex- 
tend the frequency range of radio continuum observations 
and provide important information on the structure of 
the radio continuum sources at low frequencies. In our 
610 MHz map, we only see emission from the sources 
S255-2a and S255-2b in the S255-IR cluster. These two 
sources are not resolved in the 610 MHz map, however, 
they appear as two distinct sources in our higher res- 
olution 1280 MHz map. We measured flux densities 
and sizes of the radio continuum sources using a two- 
dimensional elliptical Gaussian model via the AIPS task 
JMFIT; the results are given in Table [3j In Section |] 
we compare our data with observations at other wave- 
lengths and discuss possible relations between the radio 
continuum and molecular emission. 

Although at the VLA only spectral line observations 
were made we combined line-free channels to form con- 
tinuum images at 23.7 GHz. These continuum maps 
are not of very high quality (the usable bandwidth was 
only about 3 MHz) and so the detection level and un- 
certainties just are not that good. We detected S255-1 
and S255-2c. Their measured fluxes are 19 ± 2 mJy and 
3.3 ± 0.5 mJy, respectively. S255-2a and 2b are not de- 
tected. 
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Figure 12. Channel maps of the SiO (5-4) line emission in the S255IR area (contours) overlaid on the 1.1 mm continuum image. The 
numbers in the upper left corner indicate the channel velocity in kms -1 . The contour levels are (-3, 3, 5, 7, 10, 15, 20)x45 mjybeam -1 . 
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Figure 13. Channel maps of the SiO (5-4) line emission in the S255N area (contours) overlaid on the 1.1 mm continuum image. The 
numbers in the upper left corner indicate the channel velocity in kms -1 . The contour levels are (-3, 3, 5, 7, 10, 15, 20)x40 mJybeam - 1 . 
The dashed contours show negative features due to the missing flux. The crosses mark the positions of the SMA3 and SMA5 continuum 
clumps. 
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Figure 14. The spectrum of the SiO J = 5 — 4 emission from the 
S255N-SMA3 area averaged over the region 3' x 5 x 2'/ 5 centered at 
S255N-SMA3 (in units of brightness temperature — the conversion 
factor is approximately 2.2 K per Jy beam -1 ). 
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Figure 15. The maps of the blue-shifted (blue thick dotted con- 
tours) and red-shifted (red thick solid contours) SiO J = 5 — 4 
emission in the S255N-SMA3 area overlaid with the map of the 
1.1 mm continuum emission (green thin contours). The velocity 
ranges are indicated at the bottom of the plot. The contour levels 
for SiO are (-3, 3, 5, 7, 9)x0.75 Jy beam -1 kms~ . The contour 
parameters for continuum are the same as in Fig. ^1 The dashed 
contours show negative features due to the missing nux. 



Table 3 

Results of the GMRT radio continuum measurements. The 
fluxes, sizes (convolved) and position angles at 1280 MHz 
as well as fluxes at 610 MHz are indicated. 



Name 


#1280 


$max 


^min 


P.A. 


#610 




(mJy) 


II 


II 


(°) 


(mJy) 


S255-2a 


3.81 ±0.49 


8.1 


7.0 


46.4 


4.22 ±0.86 


S255-2b 


1.53 ±0.41 


7.7 


3.8 


55.9 


0.96 ±0.31 


S255-2c 


1.11 ±0.32 


7.7 


4.3 


58.9 


< 0.4 


S255-1 


20.25 ±0.46 


10.5 


5.0 


61.1 


10.48 ±0.78 



4. DISCUSSION 



4.1. 



General morphology and kinematics 

Our data indicate the presence of several new clumps 
in the mapped areas. Some of them are visible in 
the 1.1 mm continuum map. Most of these continuum 
clumps also show molecular emission (although the sets 
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v lsr (km/s) 

Figure 16. The spectrum of the SiO J = 5 — 4 emission from the 
S255N-SMA5 area averaged over the region 5'/ 5 x 5^5 shifted to the 
south by 4 ; /7 from the S255N-SMA5 central position (in units of 
brightness temperature — the conversion factor is approximately 
2.2 K per Jy beam -1 ). 
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Figure 17. The maps of the blue-shifted (blue thick dotted con- 
tours) and red-shifted (red thick solid contours) SiO J = 5 — 4 
emission in the S255N-SMA5 area overlaid with the map of the 
1.1 mm continuum emission (green thin contours). The velocity 
ranges are indicated at the bottom of the plot. The contour lev- 
els for SiO are (-3, 3, 5, 7) x 0.25 Jybeam -1 kms -1 . The contour 
parameters for continuum are the same as in Fig. ffl The dashed 
contours show negative features due to the missing nux. 

of emitting lines are somewhat different for each one). A 
comparison with our re cent NIR J. H an d K band ob- 
servations of this area ( |Ojha et al. 2011) shows almost 
no association of the newly detected continuum clumps 
with the NIR sources (Fig. [I]). This implies that these 
clumps are at a very early stage of evolution. 

A pronounced example is the new clump which we des- 
ignate as SMA4 in the S255IR area. It shows fairly 
strong continuum emission (Table [5| and is associated 
with emission in many molecular lines including N2H + , 
DCO+, DCN, H 2 CO, H 2 CS, and CH 3 OH. Interestingly, 
emission peaks of different species are located at dif- 
ferent distances from the continuum source, implying a 
chemical abundance gradient in approximately the east- 
northeast direction (Fig. |24|). There is also a moderate 
velocity gradient in this airection, that can be seen in 
the N2H + line channel maps (Fig. [5| and is also illus- 
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Figure 18. Maps of the CH3OH integrated line emission in the S255IR area (contours) overlaid on the 1.1 mm continuum image. The 
velocity ranges are indicated at the bottom of the plots. The contour levels are (-3, 3, 5, 7, 9, 12, 15, 20) x 0.3 Jybeam -1 kms -1 for the 
5i - 4 2 E, 61 - 7 2 A", 10 2 - 9 3 A+ and 14 4 - 15 3 E lines, and (-3, 3, 5, 7, 9, 12, 15, 20)x0.6 Jybeam" 1 kms" 1 for the other lines. The 
dashed contours show negative features due to the missing flux. 
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Figure 19. Maps of the CH3OH integrated line emission in the S255N area (contours) overlaid on the 1.1 mm continuum image. The 
velocity ranges are indicated at the bottom of the plots. The contour levels are (-3, 3, 5, 7, 9, 12, 15, 20)x0.6 Jybeam -1 kms -1 . The 
dashed contours show negative features due to the missing flux. 

ammonia lines probably have two v elocity components 
in S255IR-SMA1 and S255N-SMA1. [Wang et aL| pQlT ) 
also mentioned that the C 18 spectrum toward S255IR- 
SMA1 shows two peaks. The presence of two velocity 
components is also evident in the spectra shown in Fig. [7] 
The most intense emission toward S255IR-SMA1 in var- 
ious lines including the highest excitation lines, is ob- 
served at Vlsr ~ 4 — 5 kms -1 . Another component, at 
Vlsr ~ 10 kms -1 is seen in lines of lower excitation and 
has a somewhat smaller linewidth. 

In S255N-SMA1 the data also show two kinematic com- 
ponents — at Vlsr ~ 8 and ~ 10 kms -1 . The lines 



trated by the position- velocity diagram in Fig. 25 The 
total velocity change is about 5 kms -1 . One ofthe SiO 
clumps, partly overlapping with the N 2 H + emission, also 
appears to be associated with this object, being located 
at the largest dis tance from the continuum source (~ 1") 
(Figs. |2|5|10|1 2|). At the same time the C 18 emission 



is very weak. Perhaps we see here an outflow oriented 
almost perpendicular to the line of sight. This can ex- 
plain the relatively small velocity gradient and narrow 
line widths. 

The data also show a complex kinematics of some of 
the previously known clumps. As mentioned above the 
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Figure 20. Maps of the H2CO 3o,3 — 2o,2, H2CS 81,7 — 7i,6 and SO2 181,17 — 172,16 integrated line emission in the S255IR area (contours) 
overlaid on the 1.1 mm continuum image. The velocity ranges are indicated at the bottom of the plots. The contour levels are (-3, 3, 5, 7, 
9)x0.6 Jybeam -1 kms -1 . The dashed contours show negative features due to the missing flux. 
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Figure 21. Maps of the H2CO 3o,3 — 2o,2 and H2CS 81,7 — 7i 5 6 integrated line emission in the S255N area (contours) overlaid on the 
1.1 mm continuum image. The velocity ranges are indicated at the bottom of the plots. The contour levels are (-3, 3, 5, 7, 9, 12, 15, 19, 
24)x0.6 Jybeam -1 kms -1 for H 2 CO 3 ,3 - 2 ,2 and (-3, 3, 5)x0.6 Jybeam -1 kms -1 for H 2 CS 81,7 - 7i, 6 - The dashed contours show 
negative features due to the missing flux. 



with the excitation energy E u > 40 K are observed at 
Vlsr ~ 8 kms -1 while lower excitation lines have emis- 
sion peaks in the range Vlsr ~ 8 — 10 kms -1 . 

4.2. Basic physical parameters of the clumps 
4.2.1. Kinetic temperature 



We 
monia 



derive gas kinetic temperature 
observations us ing the usual 



from our am- 
approa ch de- 



scribed for e xample in jMangum et al. (1992); |Harju 
1H Jl993| >. The NH 3 (1,1) and (%2) spec- 
tra were^ fitted using the GILDAS software pack- 
age ( http: / /www.iram.fr /IRAMFR/GILDAS] ) . Results 
of the temperature estimates for the continuum clumps 
are given in Table [IJ Typical uncertainties are 20-30%. 
As mentioned above, in S255IR-SMA1 and S255N-SMA1 
we probably see secondary velocity components. The 
signal-to-noise ratios for them are very low; we show es- 
timates of the kinetic temperature for these components 
in parenthesis with question marks. 



In S255N-NH3 we obtained the temperature of about 
13 K for the 8 kms -1 component and about 23 K for 
the 10 kms -1 component. The temperature of the lat- 
ter component is practically the same in S255N-NH3 and 
in S255N-SMA1, while the temperature of the 8 kms -1 
component sharply increases toward S255N-SMA1 from 
very low values at S255N-NH3. The temperature esti- 
mate for S255IR-N 2 H+(1) is about 50 K. 

Our data set contains several methanol transitions of 
significantly different excitation. In principle they could 
be used for temperature estimates with the traditional 
method of rotational diagrams. Unfortunately in most 
cases the constructed rotational diagrams do not allow 
any reasonable temperature fit. The reason can be a 
significant deviation from LTE conditions assumed by 
the method of rotational diagrams. As mentioned above, 
several of the observed transitions may be masing. Pre- 
liminary results of non-LTE modelling ba s ed on the ap- 
proach described in Salii & Sobolev (2006); Sutton et al. 
(2004) show a fairly good agreement with the ammonia 
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Table 4 

Velocities and physical parameters (line width, linear size, kinetic temperature derived from 
ammonia observations, assumed dust temperature, mass, virial mass and mean density) of the 
millimeter wave continuum sources. 



Name 


Vlsr 


Ay 


L 


T kin (NH 3 ) 




M 


M vir 




n 




(kms- 1 ) 


(kms" 1 ) 


(pc) 


(K) 


(K) 


(M©) 


(M ) 


(cm" 3 ) 


S255IR-SMA1 


4.4 


3.3 


0.012 


240, (80 ?) 


40 


10 


14 


2 


x 10 8 


S255IR-SMA2 


9.2 


2.0 


0.031 


44 


40 


10 


13 




10 7 


S255IR-SMA4 


7.9 


2.6 


0.048 


< 20 


20 


14 


34 


4 


x 10 6 


S255N-SMA1 


8.0 


3.5 


0.030 


29, (> 200 ?) 


30 


23 


39 


3 


x 10 7 


S255N-SMA2 


8.7 


3.0 


0.028 


20 


2 


26 


4 


x 10 6 


S255N-SMA3 


9.0 


< 2 


0.037 


29 


30 


2 


< 10 




10 6 


S255N-SMA4 


7.0 


4.2 


0.020 




20 


8 


40 


2 


x 10 7 


S255N-SMA5 


9.5 


3.9 


0.048 


24 


25 


7 


80 


2 


x 10 6 


S255N-SMA6 


8.9 


4.5 


0.039 




20 


6 


80 


4 


x 10 6 



assumed 
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Figure 22. Our GMRT map of the investigated area at 610 MHz. 
It covers the extended HII regions S255, S256 and S257. Two 
compact sources in the center are associated with the S255IR and 
S255N star forming cores. 

temperature estimates (S. Salii & A. Sobolev, private 
communication) . Thus we use only these latter estimates 
here. 

4.2.2. Masses and densities 

Here we present some simple estimates of clump 
masses, sizes and densities derived from the dust con- 
tinuum emission. Our continuum flux measurements at 
two different frequencies (~ 225 and ~ 285 GHz) imply 
a spectral index of about 2 for both S255IR-SMA1 and 
SMA2. This indicates optically thick emission, which is 
unrealistic. A comparison of our measured flux densities 
at 225 GHz for SMA 1 and SMA2 with those published 
by Wang et ah] ( |2Q11| ) shows that our values are about 2 
times higher. This could be caused by different uv ranges 
of the observations; with our lower resolution data we are 
less susceptible to missing flux. Our 285 GHz data hav e 
an angular resolution closer to that of | Wang et al. (2011 ). 
Using our 28 5 GHz flux densitie s and the 225 GHz Hux 
densities from |Wang et al.| J2011| ), we obtain a spectral in- 
dex of about 4 tor these clumps, as expected for optically 
thin dust emission. This shows that even simple mass 
estimates based on continuum flux densities are highly 
uncertain and strongly depend on the iw-coverage. 

We assume the dust temperatures (Td) close to the 
kinetic temperatures derived from the ammonia observa- 



tions. For the clumps with no ammonia temperatures we 
assume T d = 20 K. In S255IR-SMA1 and S255N-SMA1 
we see components with different temperatures. For the 
bulk of the dust emission we assume the lowest temper- 
ature. In S255IR-SMA1 the ammonia emission is very 
weak and estimates of the lowest temperature are very 
uncertain. We assume Td = 40 K in accordance with the 
value of kinetic temperature f ound in single-dish obser- 
vations ( Zinchenko et aL] 2009 ) . Then, we assume a gas- 
007 



to-dust mass ratio of 100, and adopt a du s t abso rption 
efficiency following |Ossenkopf fc Henning ( 1994 ). The 
masses derived in this way are given in Tableland range 
from 2 to 34 M . 

In this table we also indicate velocities and typical line 
widths for the clumps as well as their linear sizes and 
virial masses. The velocities and line widths in most 
cases were derived from C 18 and/or C 34 S spectra. In 
some clumps (S255IR-SMA4, S255N-SMA6) these lines 
are too weak and we used data of other lines, in particular 
N 2 H+ for S255N-SMA6 and H 2 CO for S255IR-SMA4. 
For S255IR-SMA1 and S255N-SMA1 we give velocities 
and line widths of the most intense components. 

The clumps are elongated; we indicate linear sizes as 
the geometric mean between the major and minor axes. 
The virial masses are obtained in the simplest way, as- 
sumin g uniform spherical clouds as in | Zinchenko et al. 
( |1994[ ). The same geometry is used to derive mean gas 
densities. It is worth noting that for all clumps the line 
widths greatly exceed expected thermal widths which 

hints at strong turbulence in these objects. 

For the clumps investigated by Wa ng et al.| (2011) 
our mass estimates are very close to their estimates. 
The virial masses for the major clumps in these ar- 
eas (S255IR-SMA1, S255IR-SMA2 and S255N-SMA1) 
and for S255IR-SMA4 are close to the masses derived 
from the dust continuum emission (within a factor of 
~ 2). This means that these clumps are close to be- 
ing in gravitational equilibrium. On the other hand, the 
virial masses of the weaker clumps in the S255N area in 
some cases are an order of magnitude greater than the 
masses obtained from the dust emission. In principle, 
these clumps could in truly be unbound and represent 
transient density enhancements produced by turbulence. 
However, estimates of their virial masses are rather un- 
certain, too. For example, the ammonia lines in S255IR- 
SMA5 are almost 4 times narrower than the C 18 line 
which was used for the estimate of the virial mass. If 
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Figure 23. Maps of the 1280 MHz continuum emission in the S255IR (left panel) and S255N (right panel) areas (contours) overlaid on 
the 1.1 mm images. The contour levels are (3, 3.5, 4, 4.5, 5, 6, 7, 8, 10)x0.22 mJybeam- 1 for S255IR and (3, 3.5, 4, 4.5, 5, 6, 7, 8, 10, 
12, 15, 20, 25, 30)x0.22 mJybeam -1 for S255N. The dash-dotted (red in the color version) and dotted (blue in the color version) contours 
show the red-shifted and blue-shifted high-velocity CO emission, respectively. 
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Figure 24. Maps of the 1.1 mm continuum emission and integrated line emission of several molecu lar species in the S255IR-SMA4 
area (contours) overlaid on the 1.1 mm continuum image. The velocity ranges are the same as in Figs |2|8|10[l8| for those species. The 
contour levels are the same as earlier for continuum, N2H+ and CH3OH, and 1.5 times lower for SiO, DCN and DCO+ (in order to better 
show weaker features). The cross marks the position of S255IR-SMA4. The strait line in the N2H+ panel indicates the location of the 
position- velocity slice shown in Fig. [25] The dot on this line corresponds to the zero offset in Fig. |25| 
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Figure 25. The position-velocity diagram for the N2H+ emi ssion 
in the S255IR-SMA4 area along the line indicated in Fig. [24] 

we estimate the velocity dispersion from ammonia, the 
discrepancy between the virial mass and mass derived 
from the dust emission, disappears for this clump. The 
difference in the line widths can be caused by the fact 
that the C 18 emission is much more extended. Then, 
these clumps could be even colder than we assumed. We 
can conclude that the question of the stability of these 
clumps needs a further investigation. 

Our data show clusters of star forming clumps in both 
S255IR and S255N. Using the data on clump velocities 
presented in Table [4] we can estimate virial masses of 
these clusters. These estimates give M v i r ~ 300 M 
for S255IR and M vir - 100 M for S255N which is in 
reasonable agreement with their mass estimat es based on 
single-dish observations of th e dust emission (Zinchenko 
[eTaT1[2009l |Wang et al.poTT ). 



4.3. Molecular clumps without continuum counterpart 

In some cases, even quite strong molecular emission at 
velocities of the quiescent gas has no detectable contin- 
uum counterpart. For example, the ammonia emission 
in S255N is strong and extended toward the north in 
both the (1,1) and (2,2) lines (see Fig. [3}. This emis- 
sion occurs at distances up to ~ 25" (~ 0.3 pc) from 
the continuum source. At th e same time our own data 
and the results published by Wa ng et al.| ( |2011[ ) show 
that the ammonia emission is spatially close to tne red- 
shifted high- velocity CO emission lobe and is extended 
in the same direction. Moreover, a comparison of the 
NH3 (1,1) and (2,2) maps shows that the intensity ra- 
tio 7(2, 2) /I (1, 1) increases towards this CO lobe, imply- 
ing a temperature increase. These facts hint at a phys- 
ical association between the ammonia emission and the 
high-velocity outflow, even though they do not coincide 
spatially or have the same velocity. A possible explana- 
tion is that we see in ammonia (and partly in N2H + and 
DCO + ) the cold quiescent gas interacting with the high- 
velocity outflow. This interaction heats the gas, causing 
the observed temperature gradient. The absence of emis- 
sion in other molecular lines probably means that these 
molecules are frozen onto dust grains. This assumption 
is consistent with the low kinetic temperature derived for 
this area. However, the picture is further complicated by 



the fact that we see two kinematic components in am- 
monia ( Fig. [4 [) with significantly different temperatures 
(Section [3 2.1] ). An alternative explanation could be that 
one of the ammonia clumps is the driving source of the 
high- velocity outflow seen in CO. A point in favor of this 
interpretation is that the CO emission is not obviously 
associated with the SMA1 clump, while the ammonia 
emission is. The interpretation is hampered by the fact 
that we see no blue-shifted high- velocity CO emission. 

Another interesting zone of molecular emission with- 
out associate d conti nuum emission is S255IR-N2H + (1). 
I Wang et al.| ( |2011[ ) report methanol emission in the 
8_i — 7o F) transition and suggest that it is related to 
the shock heating excited by the outflow, and possibly is 
masing. However, our data show emission from many 
other lines, including N2H + , H2CO, several methanol 
transitions, and probably NH3 and DCO + . This seems 
to rule out the hypothesis of specific excitation of the 
8_i — 7o E methanol line and indicates the presence of 
an object with rather unusual properties. It is notewor- 
thy that we detect no C 18 or C 34 S emission in this area. 
At the same time the temperature here derived from am- 
monia, is rather high (~ 50 K) which implies a presence 
of a heating source. This area apparently also needs a 
f urther investigation . 



Cyganowski et al. (2007) and Wang et al. (2011) men- 
tioned one more area of a strong molecular emission with- 
out a continuum counterpart shifted to southwest from 
S255N-SMA1 along the direction of th e outfl ow. It is 
seen also in several our maps (e.g. Figs. |ll|21| ). 



4.4. Molecular outflows 

A description of high-velocity outflows in the S255IR 
and S255N areas based on CO J = 2 — 1 obser vations 
with the SMA was given by | Wang et al.| Q2011| ). They 
identified one outflow in the S2551R area and at least one 
in S255N, mentioning a rather complicated structure for 
the high- velocity CO emission here, and not excluding 
the possibility of multiple outflows. Our data — in par- 
ticular the SiO observations — shed new light on this 
question. First, the SiO data indicate that there are 
probably outflows originating from the weaker and pre- 
sumably less-evolved clumps. As mentioned above, we 
almost certainly see outflows associated with the S255N- 
SMA3 and S255N-SMA5 clumps. An outflow related to 
the S255IR-SMA4 clump also seems probable. None of 
these objects show NIR, radio continuum, or maser emis- 
sion, which implies a very early evolutionary stage. 

Basic physical properties of outflows are usually de- 
rived from obser vations of CO and — pre f erably — 13 CO 



1992 Henning 



and C 18 (e.g. JCabrit fc Bertout||l990, 
|et al.|2000[|Zinche nko 2002|7Tn" our case we see the high- 
velocity CO J = 2 - 1 emission towards S255N-SMA3 
and S255N-SMA5. However neither 13 CO nor C 18 are 
detected in the velocity ranges of the SiO outflows. We 
also do not see an emission in the lines of other poten- 
tial outflow tracers (e.g. methanol) at these velocities. 
Therefore we use CO data to estimate physical parame- 
ters of these outflows assuming a small optical depth in 
the CO line wings. Most probably this assumption is in- 
correct, which means that our estimates of the mass and 
mechanical parameters represent lower limits. Following 
the procedures outlined in the mentioned works we de- 
rive mass, momentum, energy, size, age, mass loss rate 
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and mechanical force for these outflows (Table |5j). The 
inclination angle was assumed to be 45°. The outflows 
appear to be very young (a few hundred years). The mass 
loss rate and mechanical force are significantly lower than 
in S255IR-SMA 1 and S255N-SMA 1, and comparable to 
those in S255S ( [Wang et al.|[2QTT| ). 

We cannot derive parameters lor the possible outflow 
related to the S255IR-SMA4 clump since we cannot iden- 
tify the relevant CO emission. 

4.5. Chemical features 

One of the most remarkable features in our molecu- 
lar maps is the absence of N2H + emission from S255IR- 
SMA1. Any N2H + emission from this clump is at least 
an order of magnitude weaker than that from SMA2. 
This fact cannot be explained by excitation effects and 
most probably implies a corresponding difference in the 
N2H + column density and the N2H + abundance, taking 
into account the similar total gas column densities for 
these clumps ( |Wang et al.|201l[ ). A drop in N 2 H + abun- 
dance in luminous massive cores was seen in our earlier 



quet et al. 



single-dish observations (Pirogov et al. 2007 Zinchenko 



2011 Reiter et all 2011 Cortes 2011 ). The 



et al. 2009 ) and reported recently by otEers (e.g ., Bus- 

T 



usual explanation lor 1N2-H+ depletion is its destruction 
by CO molecules released from dust grains at rising tem- 
perature. However, the C 18 maps do not support this 
hypothesis. The C 18 intensities are similar for SMA1 



and SMA2, if not stronger towards the latter (Wang et al. 



2011 and our data) 

It is probably not coincidental that the DCO + dis- 
tribution follows that of N2H + , as mentioned above. 
Both species apparently avoid ionized regions. This is 
also true for S255N and is consistent with our previous 
suggestion that the primary destruction mechanism for 
N2H + in such objects is dissociativ e recombination due 



to an enhanced ioniza tion fraction (Pirogov et al. 2007 
Zinchenko et al.|20 09 ) . At the same time in S255JN there 
is a noticeable difference between the DCO + and N2H + 
maps: there is no DCO + emission toward the bright- 
est N2H + peak. In principle this might be caused by 
temperature distribution in this area — an enhanced 
temperature would reduce the abundance of deuterated 
molecules. It is worth noting that the NH3 (2,2) emission 
has a peak here while emission peaks of other deuterated 
species (DCN and DNC) are also shifted from the N2H + 
peak. 

Another interesting feature is the rather strong DCN 
emission from several clumps, including S255IR-SMA1, 
which is apparently hot, as mentioned above. Although 
there are observations of rather high DCN abundances in 



warm (50-7 5 K), dense (n ~ 1Q T cm 3 ) clumps ( |Leurini 



|et al.|[2006 Parise et al.||2009l) which were successfully 



explained b y chemical models ( [Roueff et al.||2007| |Parise| 
|et al.|2009j ) it is hard to expect the same at the tempera- 
ture of ~ 150 K derived for S255IR-SMA1. Recent mod- 



elling by |Albertsson et al. ( |201l[ ) shows that DCN/HCN 
abundance ratio sharply drops at temperatures > 80 K. 
Possible explanations for a high DCN abundance in this 
clump might presume a significant amount of relatively 
cold gas in this clump or a very young age, with insuf- 
ficient time to change the isotopic ratio. The required 



peak of the DCN emission is noticeably shifted from the 
continuum peak and SiO emission peak. Most probably 
DCN traces cold material adjacent to the hot core. The 
density of this material should be high in order to explain 
the observed intensity ratio of the DCN J = 3 — 2 and 
J = 4 — 3 transitions. Our temperature estimates from 
the ammonia observations show a presence of a relatively 
cold component in S255IR-SMA1. 

In S255N the DCN emission is noticeably stronger than 
in S255IR; we also detected DNC h ere. In the framework 
of time-dependent chemistry (e.g., Hatchell et al.|p~998|) 



this may indicate that S255N is younger than Sz551R. 
This is also consistent with a lower temperature of S255N 
in respect to S255IR as found from ammonia observa- 
tions. However we need estimates of the DCN/HCN 
abundance ratios. 

As mentioned above we did not detect N2D + in ei- 
ther S255IR or S255N. An upper limit (3a) on the 
N2D + brightness implies a limit on the intensity ratio 
7(N 2 D+)/7(N 2 H+) < 0.1 for the brightest N 2 H+ emis- 
sion peaks. In cold dark clouds this ratio can reach the 
values of ~ 0.25 but drops to < 0.1 at temperatures 
> 25 K (e.g. |Emprechtinger et al. 2009). In our case, 
when temperature estimates lor the l^right N 2 H + emis- 
sion peaks are available, they show higher temperatures. 
Therefore, the non-detection of N2D + is consistent with 
our temperature estimates. 

Some other chemical peculiarities related to molecular 
emission not associated with a noticeable dust continuum 
have been mentioned above (in particular the ammonia 
source S255N-NH3 in the northern part of the S255N 
area, some of the N2H + sources, etc.). 

4.6. Properties of the radio continuum sources 



time can be less than 10 4 years as d iscuss ed by Hatchell 
et al. (1998). An inspection of Figs. 8|10 shows that the 



Our GMRT (Figs. [23]) and available VLA archival data 
show the presence oiseveral centimeter and decimeter 
wave continuum source s. The sources S255-2 a and S255- 
2b identified earlier by |Snell fc Ballyl (|1986| ) are associ- 
ated with NIR objects ( |Ujha et al.||20ll| ) that are clas- 
sified as B-type stars and apparently represent compact 
H11 regions excited by these stars. It is interesting that 
we do not see any molecular or dust emission associated 
with these sources. The high- velocity CO red-shifted 
lobe is close to S255-2a but it is unclear whether there is 
a physical association. Possibly the parental cloud mate- 
rial of these stars has already been dispersed, indicating 
a relatively large age. 

The source S255-2c is associated with S255IR-SMA1. 
It has extensions in the directions of both red-shifted and 
blue-shifted lobes of the molecular outflow, more pro- 
nounced at low frequencies. This indicates a gradient in 
the emission measure. Apparently at least part of the 
emission is associ ated with ionized gas in the outflow. 
Ojha et al] ( |2011| concluded that this source is younger 
than S255-2a and 2b. This is consistent with our data, 
which indicate the presence of a hot molecular core. 

In S255N the low frequency radio continuum maps also 
show an extension in the direction of the outflow lobe, 
as in S255IR. Presumably this also is an indication of 
emission from ionized, outflowing gas. 

The spectra of all the sources in the range from 
610 MHz to 23.7 GHz are plotted in Fig. [26] using th e flux 
densities from t he p resent observ ations, data from |SneTI 
fc Bally| ( [19861 ) and |Ojha et al.| ( [20TT] ). The flux den- 
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Table 5 

Parameters of the outflows in S255N-SMA3 and S255N-SMA5 (mass, momentum, energy, size, 
age, mass loss rate and mechanical force). 



Name 


M 


P 


E 


Size 


t 


M 


F 




(M ) 


(M© kms -1 ) 


(erg) 


(pc) 


(yr) 


(Moyr- 1 ) 


(M kms" 1 yr -1 ) 


S255N-SMA3 


0.003 


0.15 


8 x 10 43 


0.009 


200 


2 x 10~ 5 


8 x 10- 4 


S255N-SMA5 


0.012 


0.36 


10 44 


0.012 


400 


3 x 10- 5 


9 x 10- 4 




1000 



Frequency (MHz) 



Figure 26. Spectra of S255-1 (triangles), S255-2a (squares), S255- 
2b (stars) and S255-2c (circles). The data at 610 MHz, 1280 MHz 
and 23. 7 GH z are from the pr esent work, at 5 GHz from [Snell &:| 
i Ballv l l| l98 6h at 1 5 GHz from |Snell Bally| fl986] ) (S255-1) and 
[Djna et al. (20lTJ). 



sity a t 15 G Hz is different from that presented by Ojha 



et al. (2011) since the primary beam correction was not 
applied in that paper. The noise of the map at 610 MHz 
is used as upper limit for the source S255-2c and the er 



rors at 15 GHz are considered to be 15% (Ojha et al. 
2011[ private communication). It should be noted that 



these measurements are obtained from radio continuum 
maps made with different approache s and with differ- 
ent angular resolutions. For example, |Ojha et aT] ( |2011[ ) 
estimated positions and flux densities by fitting ellipti- 
cal Gaussian model components, assuming the sources 
to be point-like, where as the total flux density reported 
by|Snell&Bally|([T986l) was obtained by integrating over 
the region within the 3cr limit. 

The structure of S255N is quite comple x. The high 
resolution (~ 1 arcsec) 3.6 cm VLA map by Cyganowski 



et al. (2007) resolved the UC Hn region into three com- 
ponents, with an arc, a point source, and an extended 
tail, which add uncertainty in the estimation of fluxes. 
However, the general shape of the SED in Fig. 26 is con- 
sistent with free-free emission from thermal plasma. 

The spectra show that the emission of S255-2a is op- 
tically thin throughout this frequency range, while the 
other sources become optically thick at the lowest fre- 
quencies. An expected flux for S255-2a at 23.7 GHz is 
about the same as for S255-2c. However S255-2a is more 
extended (Table |3| and therefore its brightness should 
be lower than that of S255-2c. This can be the reason 
for its non-detection at this frequency. S255-2b is notice- 
ably weaker and apparently below our detection limit at 
23.7 GHz. 

Using the formulas for free-free emission from |Lang| 
fl!999 j and assuming an electron temperature of l(rK 



we can estimate the emission measure for these objects 
For S255-2a with a turnover frequen cy of vj_ S 60 MHz 



a et al. 
derived 



(2011) 



we obtain EM < 10° pccm _t) 
from their measured fluxes and sizes 
1.5 x 10 6 pccm -6 for this source, in agreement with our 
estimate. For S255-2b, S255-2c and S255-1 the turnover 
frequency seems to be in the range of 2-3 GHz, which 
implies EM - (1 - 2.5) x 10 7 pccm" 6 . For S255-2c 
and S255-1 these estimates ar e reasonably close to those 
derived by lOjha et al. I (|2bll| ). However, in the case of 
S255-2b, |Ojha et aLI^Uirf obtained a significantly lower 
estimate lor the emission measure which corresponds to 
the turnover frequency of about 1 GHz. The reasons 
for this discrepancy are not clear. Perhaps the flux of 
S255-2b at 5 GHz is overestimated. 

5. CONCLUSIONS 

1. We detected several new clumps in the S255IR and 
S255N areas by their millimeter wave continuum emis- 
sion. They are also detected in several molecular lines. 
Temperatures of these clumps found from ammonia ob- 
servations are ~ 20 K and their masses are estimated at a 
few solar masses. These clumps have almost no associa- 
tion with NIR or radio continuum sources, which implies 
a very early stage of evolution. At the same time, our 
SiO and other molecular data indicate the presence of 
high-velocity outflows related to some of these clumps. 
The outflows associated with S255N-SMA3 and S255N- 
SMA5 are apparently very young — a few hundred years 
only. The line widths of the clumps greatly exceed the ex- 
pected thermal widths, suggesting significant turbulence 
in these objects. 

2. In some cases there is strong molecular emission at 
the velocity of the quiescent gas, yet with no detected 
continuum counterpart. The nature of these sources is 
not clear. Some of them (e.g., the ammonia source in 
the northern part of S255N) may represent very cold gas 
not detectable in the dust continuum. Emission from CO 
and other molecules would be absent in this case, because 
they would be frozen onto dust grains. This assumption 
is consistent with the low kinetic temperature derived for 
this area from the ammonia observations. 

3. N2H + and DCO + apparently avoid ionized regions. 
There is no detectable N 2 H + emission associated with 
the SMA1 clump in S255IR (which contains a hot core 
and also ionized gas) implying a significant drop in the 
N 2 H + abundance. There is no sign that this is caused 
by CO evaporation. 

4. We detected rather strong emission in lines of 
deuterated species — DCN and DCO+ — in both S255IR 
and S255 N (in S255N the DC N J = 3 - 2 line was seen 
earlier by [Cyganowski et al.||2007|). In S255N a DNC 



emission line is also detected. In S255IR the DCN emis- 
sion is also quite strong in the vicinity of the hot core. 
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This implies either the presence of a significant amount of 
cold gas or a very young age for this hot core, insufficient 
to have changed the isotopic ratio. 

The DCO + distribution is significantly different from 
that of DCN, and qualitatively resembles that of N2H + . 
In particular, both species avoid ionized regions, and 
both are observed towards the S255N SMA6 clump (in 
contrast to most other molecules). 

5. There is no molecular material associated with the 
S255-2a and S255-2b radio continuum sources, which are 
compact Hn regions excited by B type stars. Possibly 
the parental material of these stars has already been dis- 
persed, indicating a relatively large age for these objects. 
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